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Abstract 

The general and explicit relation between the phase time and the dwell time for quantum tun- 
neling of a relativistically propagating particle is investigated and quantified. In analogy with 
previously obtained non-relativistic results, it is shown that the group delay can be described in 
terms of the dwell time and a self-interference delay. Lessons concerning the phenomenology of the 
relativistic tunneling are drawn. 
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To obtain the definitive answer for tlie time spent by particle to penetrate a classically 
forbidden region delimited by a potential barrier [l, 2, [3, 4], people have tried to introduce 
quantities that have the dimension of time and can somehow be associated with the passage 
of the particle through the barrier or, strictly speaking, with the definition of the tunneling 
times 5], y, 3, E y, 11, 12 1 . Tunneling is a general feature of wave equations and 
may be very counterintuitive when compared with the evolution of propagating waves. It 
occurs when a wave impinges on a thin barrier of opaque material and some small amount 
of the wave leaks through to the other side. Even so, it has been mainly discussed for 
the Schroedinger equation due to the shocking contrast between classical and quantum 
particles. In all cases described by the non-relativistic (Schroedinger) dynamics ll|], the 
pulse (wave packet) that emerges from the tunneling process is greatly attenuated and front- 
loaded due to the filter effect (only the leading edge of the incident wave packet survives 
the tunneling process without being severally attenuated to the point that it cannot be 
detected). If one measures the speed by the peak of the pulse, it looks faster than the 
incident wave packet. Moreover, since the transmission probability depends analytically on 
the momentum component k {T = T{k)), the initial (incident wave) momentum distribution 
can be completely distorted by the presence of the barrier of potential. 

In what concerns the momentum distribution distortion and the precise computation of 
phase times, considering the relativistic tunneling dynamics in terms of the Dirac/Klein- 
Gordon wave equation allows for circumventing such difficulties. Even though in the non- 
relativistic framework lOj, a quite elegant study is performed so to overcome the above 
mentioned misunderstanding of the tunneling time definitions. Indeed some authors consider 
difficult and perhaps confusing the treatment of all interactions of plane waves or wave 



packets with a barrier potential using a relativistic wave equation [13|, [Ij, [151, |16j. This 
is because the physical content depends upon the relation between the barrier height Vq 
and the mass m of the incoming (particle) wave, beside of its total energy E. In some 
previous analysis we have demonstrated with complete mathematical accuracy that, 
in some limiting cases of the relativistic (Klein-Gordon) tunneling phenomena where the 
relativistic kinetic energy is approximately equal to the potential energy of the barrier, and 
mcL/h << 1, particles with mass m can pass through a potential barrier VJ, of width L with 
transmission probability T approximately equal to the unity (total transmission). 

Differently from other previous (non-relativistic) tunneling analysis, the original momen- 
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turn distribution is kept undistorted and there is no filter effect. The tunnehng time is 
then computed for a completely undistorted transmitted wave packet, which legitimizes any 
eventual accelerated transmission llTf. 

^.„,.g back to the first attempt of evaluating th. p.oble™. K.eiu Q eous.de.ed the 
reflection and the transmission of electrons of energy E incident on the potential step V{x) = 
Q{x)Vo ill the one-dimensional time- independent Dirac equation which can be represented 
in terms of the usual Pauli matrices [l^ byjssi 

[ctV"5, -{E- e(x)K) - cr'^m] (t){k,x) = 0, (from this point c=h=l). (1) 

which corresponds to the reduced representation of the usual Pauli-Dirac gamma matrix 
representation obtained when the spinorial character is neglected (1 + 1 dimensional Dirac 
equation). The physical essence of such a theoretical configuration lies in the prediction that 
fermions can pass through large repulsive potentials without exponential damping, in a kind 



of (Klein) tunnelling phenomenon IJ] which follows accompanied by the production of a 



particle-antiparticle pair inside the potential barrier. It is different from the usual tunneling 
effect since it occurs inside the energy zone of the Klein paradox jisl, \\\ . 

Taking the quadratic form of the 1 + 1 dimensional Dirac equation, we obtain the Klein- 
Gordon equation for the time-like component V{x) of a Lorentz four-vector potential, 

{E - V{x)f (l){k, x) = {-dl + m') (Pik, x), (2) 

which, from the mathematical point of view, due to the second-order spatial derivatives, 
has boundary conditions similar to those ones of the Schroedinger equation and leads to 
stationary wave solutions characterized by a relativistically modified dispersion relation. 

All these proposals for computing how long a particle takes to tunnel through a potential 
barrier have led to the introduction of several transit time definitions, among which, in spite 



of no general agreement [ll|, the so called phase time |20| ( grou p delay) and the dwell 
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I22II for non-relativistic 



time have an apparently well established quantified relation 
Schroedinger equation solutions. In this manuscript, we extend such results to the one- 
dimensional scattering potential configuration described by Klein-Gordon equation solutions. 

Let us then depict the three potential regions by means of a rectangular potential barrier 
V{x), V{x) = K if < a; < L, and V{x) = if x < and a; > L. Differently from 
the non-relativistic (Schroedinger) dynamics, we observe that the incident energy can be 
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divided into three zones. The above barrier energy zone, E > Vo + m, involves diffusion 
phenomena of oscillatory waves (particles). In the so called Klein zone IJ, llSl], E < Vq —m, 
we find oscillatory solutions (particles and antiparticles) in the barrier region. In this case, 
antiparticles see an electrostatic potential opposite to that seen by the particles and hence 



they will see a well potential where the particles see a barrier 
Vq ~ m < E < Vq + m, for which only evanescent waves exist 
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25|. The tunnelling zone, 



27( 1 in the barrier region, is 



that of interest in this work. By evaluating the problem for this tunneling (evanescent) zone 
assuming that (f){k, x) are stationary wave solutions of Eq. ([2]), when the peak of an incident 
(positive energy) wave packet reach the barrier x = at t = 0, we can usually write 

(j)i{k,x) = exp[i k x] + R{k, L) exp[—i k x] x<0, 
(j)2{k,x) = a{k) exp [—p{k) x] + I3{k) eyip[p{k) x] < x < L, 
(ps{k,x) = T{k,L)exp[ik{x — L)] x > L, 

where the dispersion relations are modified with respect to the usual non-relativistic ones: 



= E^ - and p{kf 



m 



To establish a correspondence with the non-relativistic (NR) solutions, it is convenient 



to define the kinematic variables in terms of the following parameters: w = \^2mVo, v = 
Vo/m = w^/2m^, and n^(fc) = k^/w'^ = Enh/Vq. The parameter w corresponds to the same 
normalizing parameter of the usual NR analysis where fc^ = 2mENB.- The above mentioned 
relation between the potential energy and the mass m of the incident particle is given 
by the parameter v. Finally, n^(fc) represents the dependence on the energy for all the 
results that will be considered here. After simple mathematical manipulations, it is easy to 
demonstrate that the tunneling zone for the above form of the Klein- Gordon equation ([2]) 
is comprised by the interval {n^ik) — vjiy < 1 for which n^(fc) might assume larger values 
{n^{k) >> 1), in opposition to the NR case where the tunneling energy zone is constrained 
by < n'^{k) < 1). We shall observe that such a peculiarity has a subtle relation with the 
possibility of superluminal transmission through the barrier. The limits for NR energies 
[k^ « m? and V «m) are given by i;?t, << 1 and v/n « 1, which, as we have indicated 
in a previous analysis [17| , reproduces the transmission and delay results of the Schroedinger 
equation. 

The stationary phase method can be successfully applied for describing the movement of 
the center of a wave packet constructed in terms of a symmetrical momentum distribution 
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g{k — ko) which has a pronounced peak around k^. By assuming that the phase that char- 
acterizes the propagation varies smoothly around the maximum oi g{k — k^), the stationary 
phase condition enables us to calculate the position of the peak of the wave packet (highest 
probability region to find the propagating particle). With regard to the standard one-way 
direction wave packet tunneling, for the set of stationary wave solutions given by Eq. ([3]), it 



is well-known 
in terms of 



29| that the transmitted amplitude T{n,L) = |T(n, L)| exp L)] is written 



\Tin,L)\ = \l+ / smh'[pin)wL]} , (3) 



where we have suppressed from the notation the dependence on k, and 

ip{n, L) = arctan < tanh [p(ri) wL\ > , (4) 

I 2np(n) J 

for which we have made explicit the dependence on the barrier length L (parameter wL) 



and we have rewritten p{k) = w p{n), with p(n)^ = y/1 + 2n'^v — {ri^ — v /2). 

The additional phase ^p{n, L) that goes with the transmitted wave is utilized for calcu- 
lating the transit time t*-'^-' of a transmitted wave packet when its peak emerges at x = L, 

dk dn{k) dip{n, L) L 1 dip{n,L) 
~ dE{k) ~dk d^ ~vwT d^ ' 

evaluated at k = k^ (the maximum of a generic symmetrical momentum distribution g{k — ko) 

that composes the incident wave packet). By introducing the c/a55^ca/ traversal time defined 

as r^-^^ = L{dk/dE{k)) = L/v, we can obtain the normalized phase time, 

(6) 



_ f{n,L) 



T{k) g{n,L)' 
where 

f{n, L) = 8n' [(2 + Sn't; + u") - {An'' + 3u) Vl + 2n^v\ 

.V,. . o ON / TT- , o Ni Sh(p(n)wL)Ch(p(n)wL) 

p[n)wL 

g{n,L) = 16n' [2 (1 + 2n'v) - Vl + 2n^v {2n' + v)] 



+ 2 [(4 + 8n'v + v'') Vl + 2n'v - Av {1 + 2n''v)] Sh{p{n)wL)\ 

with Ch{x) = cosh (x) and Sh{x) = sinh (x). 

Turning back to the main point, could one say metaphorically that the particle repre- 
sented by the positive energy incident wave packet spend a time equal to t^^ip inside the 
barrier before retracing its steps or tunneling? The answer is in the definition of the dwell 



time for the relativistic colliding configuration which we have proposed. In quantum me- 
chanics, using steady-state wave functions, the average time of residence in a region is the 
integrated density divided by the total flux in (or out) and the lifetime is defined as the 
difference between these residence times with and without interactions. In non-relativistic 
quantum mechanics, the dwell time is a measure of the time spent by a particle in the barrier 
region regardless of whether it is ultimately transmitted or reflected I2I, 

1 



dx\<p2{k,x)\"" 



(7) 



Jin Jo 

where jin = k/m is the flux of positive energy incident particles and (l)2{k, x) is the stationary 
state wave function inside the barrier. 

In terms of the redefined parameters n and v, the explicit expression for the dwell time 
normalized by r(fc) is given by 

(8) 



foin,L) 

T{k) 



9D{n, L) 



where 



1 



n 



p{ny 
2VI + 



1 + 



\ Sh{p{n)wL) Ch{p{n)wL) 



p{n) 
Sh{p{n)wLf 



p{n)wL 



An^pinY 

To derive the relation between the dwell time and the phase time, we reproduce the 
variational theorem which yields the sensitivity of the wave function to variations in energy. 
Following from the Smith derivation 3l|] for the non-relativistic Schroedinger equation, here 
we also have the eigenvalue equation 



(ido-E) (t>{k,x) = 0, 



and its first derivative with respect to 



i^^o-E)^^-<Pik,x) = 0. 



(9) 



(10) 



After some simple mathematical manipulations and the substitution of the Klein-Gordon 
equation ([2]), the second derivative can be written as 



0. 



By following the one- dimensional analysis here considered, it is easy to find that 

^^^t_^t^^l = 2(E-VUU= — (^^-6^^\ 
dEdx^ ^ dx^dE ^ "'^ ^ dxXdEdx dEdx T 



(11) 



(12) 



where we clearly notice the presence of — K, in place of m of the non-relativistic result ([7]) 



3l| . Upon integration over the barrier length we find 

^ ^ d(t) a^t 



2{E-V,) [ dx\Mk,x)\' 
Jo 



(13) 



dE dx dEdx J 

In front of the barrier (x < 0), the wave function consists of an incident and a reflected 
component given by (f)i{k,x), and behind the barrier (x < L), there is only the transmitted 
wave 03 (fc, x) (see Eq. ([3])). Under these conditions we evaluate the right-hand side of Eq. (fT3l) 

as 

dip{k,L) Im[Ry 



2ik 



+ 1 



dA: 



k 



dk_ 
dE' 



Since + ITI 



1, Eq. (]T3|) becomes 



(S-K) / dx|,/.2(A:,x)|^ 

JO 



which gives 



T{k) 



Ad) 



dLp{k,L) k 
-dE- + e'"^^"^^ 



1 Im[R] 



r(fc) 



where we have introduced the re-scaled dwell time, 



t 



E-V, 



dx\(j)2ik,x)\'^ 



(14) 



(15) 



(16) 



(17) 



m k 

which can be related to the correct definition of the probability density for the Klein-Gordon 
equation, 

jo = [ dx [4{docl)2) - {docPD'l^^]/ t^R^ = {h/kn) = (jo/k), (18) 
Jo 

and leads to the usual definition t^°^ in the non- relativistic limit h-* 1). In Eq. ([I7D, the 

squared modulus of the wave function transforms as a Lorentz scalar, E — Vq transforms as a 
time-like component, and the integrand dx as well as k transform as space-like components of 
a Lorentz four- vector. It means that tlf ■* has the correct Lorentz character since it transforms 



as a time-like component, which does not occur for t'-^^ of Eq. (I?]). 

nn 

As in the non-relativistic case |10l . l31| , the first term of Eq. (fT6ll corresponds to the phase 
time or the aforementioned group delay. The second term comes from the explicit compu- 
tation of the dwell time. However, the presence of the multiplicative factor in Eq. f|T7|) 
introduces some novel aspects in the interpretation of the additional term —Im[R]/E as 
a self-interference term which comes from the momentary overlap of incident and reflected 
waves in front of the barrier [s]. As we can observe in the example illustrated in the Fig. [H 
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and by the usual definition ([7j), the dwell time is always positive. The re-scaled dwell time 
modulated by ( '^~^° changes sign when the total energy E equalizes the potential energy 
Vq-. an energy region comprised by the tunneling energy zone of the Klein-Gordon equation. 
Consequently, differently from the results we get from the non-relativistic analysis, the dwell 
time is not obtained from a simple subtraction of the supposed self-interference delay t^-^^^^ 
from the phase time that, in some circumstances 17|, |29|] describes the exact position of 
the peak of the scattered wave packets. Such results give a complete description of the 




1.0 1.5 



Nomnlized Energy n(k) ~ k'/2mV 



2.0 2.5 3.0 

Normalized Energy n(k) ~ li/2mV 
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FIG. 1: Delay times calculated from the dynamics of the Klein-Gordon equation: Phase Time 
(Dash-dotted line), Dwell Time (Dashed black line), Self-Interference term (Solid line), and the 
re-scaled Dwell Time (Dashed red line). In fact the tunneling region is comprised by the interval 
{•n? — v/2y < 1, > 0. Here we have adopted the illustratively convenient value of wL = 27r with 
u = 5, in comparison with the non-relativistic results parameterized by u — > 0. . 



the tunneling zone, Va — m < E < Vq + for which only evanescent waves exist [26|, l27j . 
several times ignored in the analysis of relativistic tunneling. For the evanescent tunneling 
zone, the Dirac equation and its quadratic form (namely, the Klein-Gordon equation) leads 
to the same results when we apply the (evanescent) tunneling time definitions in which we 
are interested, the phase time and the dwell time. The evanescent tunneling zone does not 
intersect with the Klein paradox energy zone for which, at least theoretically, the possibility 
of creation/annihilation of fermionic pairs leads to the reinterpretation of the probability 
density currents, and thus to a novel interpretation of the tunneling phenomenon. Con- 
sequently, concerning the calculation of evanescent tunneling times, all the references to 
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fermionic (Dirac) and bosonic (Klein-Gordon) particles are valid, in the same sense that all 
the results derived from the non-relativistic Schrodinger equation are supposed to be valid 
for massive fermions and massive bosons. 

At least for the moment, the above (relativistic) results do not necessarily demand for a 
confront with the (non-relativistic) predictions derived from the opaque limit analysis which 
results in the filter effect and the superluminal tunneling. To clear up this assertion, it 
is convenient to recover the limiting configurations [n^ v/2^il) of some of our previous 
results 



171 ] for which the tunneling transmission probability ([3]) can be approximated by 

^ K >> 1 



lim \T{n,L) 

n2^i;/2=Fl 



1 + 



2^:^:4 



[1 + (mL)r^ , (19) 



from which, avoiding any kind of filter effect, we recover the probability of complete tunneling 
transmission when rriL << 1, once we have \T{n,L)\ ^ 1. For the correspondent values of 

the phase times we obtain |l7t |. 

ti-fi) 4 1 

lim = -oTXT^' ri'^v/2Tl, n^ u > 0, (20) 

n2-^D/2Ti r(fc) 3 1 ± 2n2 

that does not depend on rriL, and we notice that its asymptotic (ultrarelativistic) limit always 
converges to 0. Curiously, in the lower limit of the tunneling energy zone, v/2 — 1, it 

is always negative. Since the result of Eq. (!20|) is exact, and we have accurately introduced 
the possibility of obtaining total transmission {transparent barrier), our result ratifies the 
possibility of accelerated transmission (positive time values), and consequently superluminal 
tunneling (negative time values), for relativistic particles when mL is sufficiently smaller than 
1 (^ T ^ 1). By observing that the barrier height has to be chosen such that one remains in 
the tunneling regime, it is notorious that the transmission probability depends only weakly 
on the barrier height, approaching the perfect transparency for very high barriers, in stark 
contrast to the conventional, non-relativistic tunneling where T{n,L) exponentially decays 
with the increasing V^. Obviously, the above results correspond to a theoretical prediction, 
in certain sense, not so far from the experimental realization. The above condition should 
be naturally expected since we are simply assuming that the Compton wavelength {h/{mc)) 
is much larger than the length L of the potential barrier that, in this case, becomes invisible 
for the tunneling particle. In general terms, the relativistic quantum mechanics establishes 
that if a wave packet is spread out over a distance d » 1/m, the contribution of momenta 



IpI ~ m >> is heavily suppressed, and the negative energy components of the wave 
packet solution are negligible; the one-particle theory is then consistent. If we want to localize 
the wave packet in a region of space (wave packet width d) smaller than or of the same size 
as the Compton wavelenght, that is d < 1/m, the negative energy solutions (antiparticle 
states) start to play an appreciable role. The condition d < L < l/m (where d < L is 
not mandatory) imposed over a positive energy component of the incident wave packet in 
the relativistic tunneling configuration excite the negative energy modes (antiparticles) and, 
qualitatively, report us to the Klein paradox and the creation of particle-antiparticles pairs 
during the scattering process which might create the intrinsic (polarization) mechanisms for 
accelerated and/or non-causal particle teletransportation. 

To conclude, in analogy with previous non- relativistic results [l^, we have shown that 
the group delay can be described in terms of the dwell time and a self-interference delay. 
The general and explicit relation between phase times and the dwell times for quantum 
tunneling of a relativistically propagating particle was investigated and quantified. Our 
analysis corroborates with the statement of conditions for the occurrence of accelerated 
tunneling transmission probabilities at nanoscopic scale in confront with the problematic 
superluminal interpretation originated from the study based on non-relativistic dynamics of 
tunneling 17|. By eliminating the filter effect, the transmission probabilities approximates 
the unitary modulus (complete transmission through a transparent medium) . In this case, we 
have noticed the possibility of accelerated {t^ < Tj^f^j ), and eventually superluminal (negative 
tunneling delays, t^ < 0) transmissions without recurrin g to the usual analysis of the opaque 
limit (p{n)wL oo) which leads to the Hartman effect [23|. 
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